Depinning and creep motion in glass states of flux lines 
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Using dynamical computer simulation we have investigated vortex matters in glass states. A 
genuine continuous depinning transition is observed at zero temperature, which also governs the 
low-temperature creep motion. With the notion of scaling, we evaluate in high accuracy critical 
exponents and scaling functions; we observe a non-Arrhenius creep motion for weak collective pin- 
ning where Bragg glass (BrG) is stabilized at equilibrium, while for strong pinning the well-known 
Arrhenius law is recovered. In both cases, a sharp crossover takes place between depinning and creep 
at low temperatures. The possible relation between the present results and a recent experimental 
observation of a second-order like phase boundary inside the BrG phase is discussed. 

PACS numbers: 



Introduction - Depending on the strength of the random 
pinning force, the equilibrium state of the flux lines at 
low enough temperature can be either Bragg glass (BrG) 
[3, [H, where quasi long-range lattice order still survives, 
or vortex glass (VG) [3|] as random as the liquid phase. 
The competition between the vortex repulsion and the 
random pin potential builds up a highly nontrivial en- 
ergy landscape, which manifests itself drastically in dy- 
namics, e.g. vortex motion under current driving. Since 
the proposal of collective pinning theory of vortex lines 
by Larkin and Ovchinnikov Q, theoretical understand- 
ing for the nonlinear dynamics response has been ad- 
vanced [l|, 0, [H, 0, HI- The functional renormalization 
group (FRG) has also been formulated 0. For recent 
review articles see [HI 03 • 

The current-driven flux lines compose a unique system 
of D = 3 dimensions in the internal space, equal to the 
dimension d — 3 of the space where the system is embed- 
ded, and N = 2 components of displacement vector. A 
full FRG treatment is still not available for N > 1. We 
try to address the issue by computer simulations, with 
the hope that useful insights complimentary to the previ- 
ous works can be provided. It is recalled that computer 
simulations for a domain wall in a plane with N = I, 
D = 1 and d = 2 were reported recently [HI]. 

On the other hand, a new experimental finding of a 
second-order like phase boundary in H — T phase di- 
agram was reported very recently, which intersects the 
first-order phase boundary associated with the melting 
transition of the BrG [Til ]. It certainly renews interests 
on corresponding variations in dynamical properties in- 
side the BrG phase; the new phase transition was dis- 
cussed in terms of replica-symmetry breaking [15fl , which 
is perhaps best captured by dynamical responses. 

The main results of the present work are summarized 
as follows: Based on the simulation results, we have de- 
rived a scaling relation among the velocity, force and 
temperature with two universal exponents for vortex mo- 
tions around the zero-temperature depinning force. From 
the exponents and the scaling curve, the Arrhenius law 
for the creep motion with a linearly suppressed energy 



barrier appears for strong pinning strength. A non- 
Arrhenius type creep motion is observed at weak collec- 
tive pinning for which the equilibrium state is a BrG. A 
sharp crossover between the depinning and creep motion 
is also observed. 

Model and simulation details - We consider a supercon- 
ductor of layered structure with magnetic field perpen- 
dicular to the layers. The model system is a stack of su- 
perconducting planes of thickness d with interlayer space 
s. Each plane contains N v pancake vortices (PV) and N p 
quenched pins. The overdampe d eq uation of motion of 
the ith pancake at position r, is [lrj, llH, 



VP I 



' th- 



Here r\ is the viscosity coefficient. The vortex-vortex 
interaction contains two parts: an intraplane PV-PV 
pairwisc repulsion given by the modified Bessel func- 
tion U (pij,Zij — 0) = deoKo(pij I X a b), and an inter- 
plane attraction between PVs in adjacent layers given by 
U vv {fHj,zt3 =s) = ( Seo /7r)[l + ln(A Qb / S )][( Pu /2r g ) 2 -l] 
for pij < 2r g and U vv (p i j 1 z lJ = s) = (se /7r)[l + 
ln(X ab / s)][pij/r g - 2] otherwise, where e = </%/ '2ir p X 2 ab 
with X a b the in-plane magnetic penetration depth, r g — 
J^ab with £ Q b the in-plane coherence length and 7 the 
anisotropy (T5 , [l9j . Here, p^ is the in-plane compo- 
nent of position vector ry between ith and jth PVs 
and Zij is that along layer normal. The pinning po- 



tential is given by U (p 



Ji Pl -= -aA p exp[-(p ip / R p ) 2 ], 
where A p = (e d/4)ln[l + (i^/2^ b )] and a is the di- 
mensionless pinning strength. Finally, is the uni- 
form Lorentz force, and F t h is the thermal noise force 
with zero mean and a correlator (Ffi(z,t)FR(z' ,t')) = 
2T]TS afi 6(z - z')8{t - t') with a,P = x, y. We consider a 
material similar to Bi2Sr2CaCu20s with k = X ao /!;ab = 
90, 7 = 100, and d = 2.83 x 10- 3 A afc , s = 8.33 x 10- 3 A afc , 
with R p = 0.22A a fc [19]. Below the units for length, en- 
ergy, termperature, force and time are taken as X a b, deo, 
deo/ks, deo/ X a b, and rjX 2 ab /d€Q. Pinning strength is sup- 
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FIG. 1: v — F characteristics for a — 0.2. Data for T — 
0.003 obtained for L xy = 40A a t, are also included (light-blue 
squares). Inset: data for zero temperature. Error bars are 
smaller than the size of symbols. 



posed to be uniform. Periodic boundary conditions are 
put in all the three directions. The results shown be- 
low are for N v = 180, N p = 900, N z = 20, L xy = 30. 
The magnetic field is roughly B ~ 100G if we take 
X a b = 2000A [2l| . The equation is integrated by the 2nd 
order Rugen-Kutter algorithm with At = 0.002 ~ 0.01. 
All data presented below are the average over 10 sam- 
ples with different randomly distributed pins. Fixing 
the magnetic field, we have performed simulations for 
L xy = 20 and 40, and make sure that finite-size effects 
are negligible to the main, universal results. For pinning 
strength, we choose a = 0.2 and a = 0.05, for which the 
maximal curvatures of the pinning potential are 4.87 and 
1.22, above and below the tilt modulus of the flux lines 
C ~ 1.47; they are therefore expected to fall into the 
strong and weak collective pinning regimes respectively 
[23,123. 

Strong pinning - Let us begin with the strong pinning 
case of a = 0.2. The equilibrium state is a pinned solid 
(VG) as random as liquid seen from the structure factor. 
The v — F characteristics at low temperatures are de- 
picted in Fig. [T](T = in inset). A continuous depinning 
transition is observed at T — with a unique depinning 
force [24|, which can be described by v ~ A(F/F c o — 1)^ 
with F c0 ~ 0.231 ± 0.002 and (3 ~ 0.74 ± 0.02. For 
< T < 0.0005, upward-convex v — F characteristics are 
observed down to F c o; below F c q there is an extremely 
small tail which is hard to see in the present scale. 

The sharp depinning transition is clearly rounded by 
finite temperatures. In order to explore the critical prop- 
erties of the v — F characteristics at finite temperatures, 
we postulate the following scaling ansatz [2~3, [26|, [2?| 



v(T,F)=T 1 / s S{T~ 1 / 05 f) (2) 
with / = 1 — F c0 /F [28[ and the scaling function S(x) 
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FIG. 2: Scaling plot for the data in Fig[T] Inset: Temperature 
dependence of velocity for several forces around F c a- 



regular at x — and S{x) — > x' 3 as x — > +00. 

The critical force can be determined using the prop- 
erty v£T, F = F c0 ) = S , (0)T 1 /' 5 implied in Eq. © (see 
also [2l|). As shown in the inset of Fig. 2, we evalu- 
ate F c o = 0.2315 ± 0.0013, and meanwhile from the slope 
1/8 = 0.754±0.010. We then perform the scaling plot us- 
ing the scaled variables vT- 1 ' 5 and (1 - Fco/F)^ 1 ^ 6 ; 
the best collapsing of data to a single scaling curve is 
achieved when (3 ~ 1/8 — 0.754, thus determining the 
exponent (3. The values of F c q and (3 estimated from 
data at finite temperatures via the scaling analysis are 
consistent with those derived from T = 0, which can be 
taken as an evidence for the existence of scaling. Fit- 
ting the scaling curve, we obtain S{x) = O.lAx 13 + 2.6 for 
x > and S(x) = 2.6 exp(0.027x) for x < 0; the former 
covers the continuous depinning transition at T — 0; the 
latter, combined with the relation f38 = 1, indicates that 
the motion at low temperatures and forces below F c q is 
well described by the Arrhenius law, and that the energy 
barrier disappears linearly when the force is ramped up 
to F c q; the bare energy barrier is U c — 0.027. 

The same exponents and similar scaling behaviors are 
available for a = 0.4; therefore, the above properties are 
universal for strong pinning case. The creep law derived 
above confirms the a priori assumption in the Anderson- 
Kim theory [3(|, and coincides with the FRG results in 
Ref. [3l[ for a domain wall (N =1). Our results are also 
consistent with Ref. [2!| provided 6 = 1. 

Deviations from the scaling curve are observed for (a) 
F c0 /2 < F < F c0 at T > 0.015 ~ U c /2, due to extra 
thermal deformation of flux-line lattice; (b) F < F c q/2 
which may be governed by the zero-force limit; and (c) 
F ^> F c o the flux-flow regime. 

Weak collective pinning - We have performed the same 
simulations for a = 0.05, which falls into the weak collec- 
tive pinning regime. The equilibrium state is a BrG, with 
the melting temperature T m ~ 0.077 above which the 
quasi long-range order is suppressed by thermal fluctua- 
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FIG. 3: v — F characteristics for a — 0.05. Data for T — 
0.0006 obtained for L xy — 40A a 6 are also included (light-blue 
squares) . 
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FIG. 4: Scaling plot for data in Fig. |3] 



tions. Looking at the velocity-velocity correlation func- 
tion, we find in this case that the moving system shows 
good t emp oral and spatial orders characterized by Bragg 
peaks [12 HH (No moving smectic [H| was observed for 
the present parameters.) The v — F characteristics are 
shown in Fig. [31 where similar to the strong pinning case, 
a sharp crossover between the depinning and creep mo- 
tion takes place at F c0 for < T < 0.00008. The scaling 
plot is depicted in Fig. [4] The exponents are estimated 
as S = 2.3 ± 0.1 and (3 = 0.65 ± 0.01, which are different 
from those for strong pinning case. The product of the 
two exponents (38 = 3/2 deviates from unity, indicating a 
nonlinear scaling relation between temperature and force 
deviation from F c q. 

Nonlinear scaling relations have been found in CDW 
systems [26j. This behavior can be captured by an effec- 
tive potential of linear and cubic terms of displacement, 
with a small energy barrier [H,[26|]. The linear scaling re- 
lation observed for strong pinning force corresponds to an 
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FIG. 5: Vortex motions at T = 0. (a) and (b): Trajectories 
of 7 nearest-neighboring vortices on one layer for a — 0.2 at 
F = 0.24 and a = 0.05 at F = 0.027 respectively, (c) and 
(d): Trajectories of the two flux lines indicated in (a) and (b). 
The total times are 330 for a = 0.2, and 1200 for a = 0.05, 
such that flux lines in both systems move roughly the same 
distance 4A ~ 1.7ao in average. Open circles in (a) and (b) 
represent the positions of pins, while the filled circles in (c) 
and (d) represent pins that pin the flux lines. 



effective potential of linear, quadratic and quartic terms 
of displacement, with a large energy barrier 35]. It is 
observed in our simulations that the system with weak 
bare pinning strength experiences smaller energy barri- 
ers compared with that of strong bare pinning strength, 
even at the same relative force-deviation from the criti- 
cal values; an analytic derivation of the effective energy 
barrier, however, is not an easy task [T^]. While a full 
picture remains to be developed, we notice that, first, 
vortex motions mimic nucleation processes in first-order 
phase transitions: the weak/strong pinning case corre- 
sponds to a system locates at the spinodal/coexistence 
curve [351 ] ; secondly, under weak and strong pinning vor- 
tices behave similarly to CDW [26[ and domain wall [111 
respectively, due to the different ranges of correlation. 

The scaling curve in Fig. Q] is fitted well by S(x) = 
0.15 exp(0. 0345a;), which, with the nonlinear scaling vari- 
able, indicates a non-Arrhcnius creep motion for the weak 
pinning case. To the best of our knowledge, this behav- 
ior has not been reported so far. On a phenomenological 
level, it can be ca ptu red with an appropriate correlator 
of energy barriers [361 ] . 

"Microscopic " vortex motion- We have examined the mi- 
croscopic motions of flux lines. As detailed in Fig. [5] for 
forces slightly above F c q and at zero temperature (similar 
results have been obtained for F < F c q at finite temper- 
atures), flux lines in the weak pinning case move homo- 
geneously in an intermediate time scale, which guaran- 
tees the moving BrG (mBrG) order [33| ; flux lines under 
strong pinning move in an inhomogenous way, and thus 
additional dislocations are induced during the motion. 
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Even for the case of mBrG, velocity fluctuates for time 
being and between different vortices, as can be seen in 
Fig. [5ji. Namely in a shorter time scale, the motion of 
vortices in weak pinning case is also random: Flux lines 
exhibit an intermittent pattern of motion, i.e. a portion 
of flux lines move while others are almost motionless at a 
given time window, and in the next time window, similar 
situation occurs with the moving vortices found in differ- 
ent regions. The motion of the total system is thus an 
accumulation of individual, random movements, in which 
creep events dominate and the random pinning plays an 
essential role. As shown in Figs. and d the shape of 
vortex lines upon depinning depends on pin strength [23| . 

Disscussions- The strong pinning a = 0.2 and weak col- 
lective pinning a — 0.05 according to the Labusch cri- 
terion [22|, [H| result in vortex glass and Bragg glass re- 
spectively. But generally the Labusch criterion does not 
coincide with the phase boundary of Bragg glass. Further 
study is expected to resolve this point in detail. 

As can been seen in Figs. [5^ and b, during the creep 
motion vortices displace frequently into the direction 
transverse to the force and the averaged velocity. In this 
way a bypath of lower energy can be found which assists 
vortices to avoid otherwise high energy barriers. This 



may affect the averaged velocity significantly since the 
time spent in large energy barriers is huge, and is to be 
taken into account carefully in theoretical treatments. 

Beside the genuine, zero-temperature depinning tran- 
sition, we observe a quite steep crossover between the 
depinning and creep motions at low temperatures (T < 
0.0005 for a = 0.2 and T < 0.00008 for a = 0.05). The 
velocity tail below the crossover is very small as such it 
looked like a phase transition in simulations (and per- 
haps also in experiments). The crossover, however, loses 
its sharpness when temperature increases but still far be- 
low the melting temperature, which might be not quit e 
consistent with experimental phase boundary H g (T) [l4|. 
The possibility of a transition (or sharp crossover) be- 
tween the depining and flux flow regimes remains to be 
investigated within the present formalism, taking into ac- 
count the shaking magnetic field. 
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